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SUMMARY 


The  tumbling  motion  of  aerodynamically  stable  bodies  entering  planetary 
atmospheres  is  analyzed  considering  that  the  tumbling,  its  arrest,  and  the  sub¬ 
sequent  oscillatory  motion  are  governed  by  the  equation  for  the  fifth  Painleve 
transcendent.  Results  based  on  the  asymptotic  behavior  of  the  transcendent  are 
applied  to  study  (l)  the  oscillatory  behavior  of  planetary, probe  vehicles  in 
relation  to  aerodynamic  heating  and  loads  and  (2)  the  dynamic  behavior  of  the 
Australian  tektites  on  entering  the  Earth's  atmosphere,  under  the  hypothesis 
that  their  origin  was  the  Moon. 


INTRODUCTION 


When  a  body  in  space  is  separated  from  a  parent  body ,  in  most  cases  the 
separation  process  will  leave  the  body  with  a  measure  of  angular  momentum. 

The  body  then  will  tumble  at  a  constant  rate  about  its  center  of  gravity  for 
the  remainder  of  its  sojourn  in  space.  As  it  enters  a  planetary  atmosphere, 
however,  its  tumbling  rate  will  begin  to  be  affected  by  the  aerodynamic  forces 
and  moments  that  come  into  play  as  the  dynamic  pressure  builds  up.  If  the 
body's  shape  is  such  as  to  provide  a  measure  of  aerodynamic  stability,  the 
increasing  aerodynamic  moment  eventually  will  bring  the  tumbling  rate  to  zero. 
Subsequently,  the  body  will  undergo  an  oscillatory  motion  of  pqssibly  large 
but  diminishing  amplitude.  It  is  of  interest  to  analyze  the  sequence  both  in 
connection  with  the  design  of  vehicles  destined  to  make  uncontrolled  entries 
into  planetary  atmospheres  and  in  connection  with  the  study  of  the  motions  of 
extraterrestrial  objects  that  have  found  their  way  into  the  Earth's 
atmosphere. 

The  equations  governing  a  tumbling  entry  are  inherently  nonlinear  and 
hence  difficult  to  treat  analytically.  The  first  investigation  in  this  field 
therefore  was  a  numerical  study  (ref.  1)  of  a  particular  vehicle  entering  a 
specific  atmosphere  (Mars').  Subsequently,  an  analytical  treatment  of  the 
problem  was  reported  in  reference  2.  Results  of  that  study  showed  that  sim¬ 
plifying  approximations  could  be  made,  enabling  the  tumbling  motion,  its 
arrest,  and  the  subsequent  oscillatory  motion  to  be  governed  by 'a  single  dif¬ 
ferential  equation.  This  equation  was  identified  as  the  equation  for  the 
fifth  Painleve  transcendent.  Following  that  study,  another  analysis  was 


carried  out  (ref.  3)  in  which  the  functional  relationships  existing  between 
the  motion  of  the  vehicle  and  the  vehicle  and  planetary  properties  were 
deduced  from  the  asymptotic  behavior  of  the  transcendent. 

The  purpose  of  the  present  work  is  to  combine  the  analyses  of 
references  2  and  3  in  a  unified  treatment  and  to  apply  the  results  to  two  dif¬ 
ferent  types  of  problems.  In  the  first,  the  tumbling  entry  of  a  probe  vehi¬ 
cle  is  considered,  and  it  is  shown  how  the  results  may  be  used  to  make  rapid 
estimates  of  the  amplitudes  of  the  oscillatory  motion  in  relation  to  aerody¬ 
namic  heating  and  loads.  In  the  second,  attention  is  directed  to  an  aspect  of 
the  research  on  the  origin  of  tektites.  On  the  basis  of  the  work  of  Chapman 
and  Larson  (ref.  4)  the  Australian  tektites  are  taken  to  be  of  lunar  origin 
and  the  consequent  dynamic  behavior  of  the  tektites  on  entering  the  Earth1 s 
atmosphere  is  examined  within  the  context  of  the  present  analysis. 


SYMBOLS 


A 

reference  area 

B 

parameter  defined 

by  equation 

(43) 

C 

parameter  defined 

by  equations 

(27) 

CD 

drag  coefficient, 

drag 

qA 

lift 

lift  coefficient, 

qA 

Cm  pitching -moment  coefficient, 

f  dependent  variable,  tan  e 

g  acceleration  due  to  gravity 

G  parameter  defined  by  equations  (27) 

Hs  stagnation -point  heat  transferred  per  unit  area 

I  pitching  moment  of  inertia  about  center  of  gravity 

Iz  moment  of  inertia  about  axis  of  spin 

J0(x)  Bessel  function  of  first  kind  of  zero  order 

Jx(x)  Bessel  function  of  first  kind  of  first  order 

l  reference  length  for  moment  coefficient  evaluation 

m  body  mass 
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integer,  denoting  number  of  complete  tumbles 

dynamic  pressure,  ^  pV2 

distance  from  center  of  planet  to  body 

dynamic -pressure  parameter,  PVj_  sin  7  ± 

time 

horizontal  component  of  flight  velocity  (sketch  (a)) 
vertical  component  of  flight  velocity  (sketch  (a)) 
flight  velocity  (sketch  (a)) 
independent  variable  (eqs.  (10)) 

axes  fixed  in  space  with  origin  at  planet  center  (sketch  (a)) 
altitude 

Bessel  function  of  second  kind  of  zero  order 
Bessel  function  of  second  kind  of  first  order 
aJQ(x)  +  bY0(x) 

angle  of  attack  in  planar  motion  (sketch  (a)) 
density  parameter  (eq.  (7)) 

flight -path  angle,  positive  when  depressed  from  local  horizontal 
(sketch  (a)) 

dependent  variable  (eq.  (l6)) 

initial  value  of  e  which  causes  body  eventually  to  trim  in  backward 
attitude 

increment  in  initial  value  of  e  measured  from  ec  (eq.  (19)) 
angle  of  pitch  measured  from  axis  fixed  in  space  (sketch  (a)) 
angle  of  pitch  measured  from  local  horizontal  (sketch  (a)) 

initial  value  of  x  (eqs.  (10)),  J-  J  q^  Cm 

ablation  parameter  (eq.  (63)) 
atmospheric  density 
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CJ 
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atmospheric  density  at  surface  of  planet 

resultant  angle  of  attack  in  nonplanar  motion  (sketch  (h)) 

displacement  of  body  from  fixed  space  axis,  0-0  (sketch(a)) 
velocity  of  rotating  liquid  (sketch  (e)) 


value 

value  of  quantity  at  first  peak  of  oscillatory  motion 

maximum  value  of  quantity 

minimum  value  of  quantity 

envelope  of  oscillatory  function 

quantity  evaluated  at  maximum  heating 


angular 


angular 

_d_ 

dt 

_d 

dx 

initial 


(  ) 
(  ) 


ANALYSIS 


The  analysis  leading  to  a  differential  equation  that  can  be  said  to 
characterize  the  tumbling  entry  problem  was  initially  presented  in  reference  2. 
That  analysis  will  be  given  here  again  in  a  slightly  amended  form  to  bring  out 
its  relevance  to  the  study  of  the  motions  of  extraterrestrial  bodies  as  well 
as  of  vehicles.  As  before,  in  order  to  reduce  the  problem  to  manageable  pro¬ 
portions,  it  is  assumed  at  the  outset  that  (l)  the  rotation  of  the  planet  and 
of  its  atmosphere  may  be  neglected;  (2)  the  motion  is  planar;  and  (3)  the 
acceleration  due  to  gravity  is  constant.  Further  assumptions  and  approxima¬ 
tions  will  be  introduced  as  necessary. 


Equations  of  Motion 


Under  the  above  assumptions ,  the  equations  governing  the  body*s  path  and 
its  motions  about  that  path  may  be  written  as 


-mV  -  CpqA  +  mg  sin  7=0 
mVj  +  CLqA+  m[(v2/r)  -  g]  cos  7  =  0 

10  -  qAZCm  =  0 


(1) 


k 


The  angles  a,  7 ,  £),$,©  are  defined 
in  sketch  (a) . 

Simplified  angle -of -attack 
equation.-  Numerical  solutions  of  equa¬ 
tions  (1)  indicate  that  after  the  body 
enters  the  planetary  atmosphere,  there 
is  always  an  interval  over  which  the 
flight -path  angle  y  and  the  flight 
speed  V  do  not  change  significantly. 
Since  this  is  the  interval  over  which 
any  tumbling  motion  would  occur,  it  is 
an  appropriate  approximation  to  take 


(2) 


7  «  const  =s  7j_ 

V  «  const  =  Vj_ 

Then,  since  0  =  u^r,  0  "will  be  essentially  zero,  so  that 


Sketch  (a) 


©  »  a 


(3) 


Next,  it  will  be  assumed  that  aerodynamic  damping -moment  terms  are  negligible 
over  the  range  of  interest,  so  that  the  aerodynamic  moment  in  equations  (1)  is 
a  function  of  angle  of  attack  only.  This,  in  conjunction  with  equation  (3), 
enables  one  to  consider  the  last  of  equations  (l)  independently  of  the  other 
two.  That  equation  becomes 

a  -  f-  q(t)Cm(a)  =  0  W 

Aerodynamic  restoring -moment  coefficient. -  As  a  body  tumbles ,  it  sweeps 
through  the  entire  angle -of -attack  range.  The  aerodynamic  restoring -moment 
coefficient  therefore  must  be  specified  over  that  entire  range.  Consider 
first  a  vehicle  intended  to  act  as  a  planetary  probe.  It  is  reasonable  to 
anticipate  that  heating  and  stability  considerations  will  dictate  the  choice 
of  its  shape,  hence,  that  it  will  be  short  with  a  conic  profile.  Also,  in 
order  to  minimize  the  amount  of  heat  protection  required,  it  is  advisable  that 
the  vehicle  be  statically  stable  in  one  trim  position  only.  In  reference  1  it 
was  found  that  the  latter  requirement  could  be  fulfilled  by  the  addition  of  a 
convex  afterbody  to  the  original  conic  forebody.  Inspection  of  the  experi¬ 
mental  results  collected  in  reference  1  for  a  vehicle  of  this  shape  reveals 
that  the  aerodynamic  restoring -moment  coefficient  as  a  function  of  angle  of 
attack  is  approximately  a  sine  wave.  Accordingly,  it  will  be  assumed  that  for 
the  probe  vehicle,  Cm(a)  in  equation  (k)  may  be  approximated  by 

Cm(°0  =  Cmmax  sin  a  (5) 

where  is  presumed  to  be  available  either  from  experimental  data  or, 

for  example,  from  Newtonian  impact  theory.  The  choice  of  equation  (5)  as  an 
appropriate  form  for  the  aerodynamic  restoring -moment  coefficient  may  be 
arrived  at  from  another  point  of  view.  With  the  addition  of  a  convex  after¬ 
body,  the  shape  whose  aerodynamic  restoring -moment  coefficient  is  reasonably 
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well  approximated  by  equation  (5)  ma y,  with  a  little  imagination,  be  thought 
of  as  a  modification  to  a  sphere.  Now  it  is  easy  to  see  that  a  sphere  with 
its  center  of  gravity  displaced  from  its  center  of  volume  has,  without 
approximation,  precisely  equation  (5)  as  the  form  of  its  aerodynamic 
restoring -moment  coefficient.  In  this  light,  therefore,  it  is  not  surprising 
that  equation  (5)  should  appear  as  the  appropriate  choice.  In  this  light 
also,  consider  the  tektites  found  in. southeast  Australia.  As  pointed  out  in 
reference  k,  the  primary  shapes  of  about  80  percent  ..of  them,  that  is,  their 
shapes  on  entering  the  atmosphere,  were  undoubtedly  spheres  or  spheroids. 

The  aerodynamic  restoring -moment  coefficient  for  a  perfect  sphere  is,  of 
course,  identically  zero;  but  with  a  small  amount  of  oblateness  which  may  be 
considered  as  an  equivalent  displacement  of  the  mass  center  from  the  center 
of  volume  of  a  perfect  sphere,  the  tektites  then  would  have  had  aerodynamic 
restoring -moment  coefficients  of  the  form  of  equation  (5)-  It  is  this  fact 
that  places  the  study  of  their  motions  within  the  framework  of  the  present 
analysis . 

Dynamic -pressure  history.-  Consistent  with  the  approximations  underlying 
equations  (2),  the  altitude history  of  the  body  as  a  function  of  time  is 

y  -  y±  =  -v.jt  (6) 


where 


v±  =  vi  sin  y± 

The  assumption  that  the  planet's  atmospheric  density  varies  exponentially 
with  altitude 


P  = 


poe 


-Py 


then  gives  for  the  dynamic  pressure  over  the  range  of  interest 


q(t) 


with 


(7) 

(8) 


q.  ~  —  p  y.2e“Pyi 

^-l  2  po  l 


Alternatively,  if  a  precise  time  history  of  the  dynamic  pressure  is  available, 
a  more  accurate  estimate  of  q.  and  Pv-^  may  be  obtained  by  fitting  the  best- 
straight  line  to  the  initial  portion  of  the  dynamic -pressure  history  plotted 
on  semilogarithmic  paper.  In  this  regard,  it  should  be  clear  that  equa¬ 
tion  (8)  reveals  one  of  the  apparently  more  severe  limitations  of  the  present 
analysis,  namely,  that  it  can  be  expected  to  apply  only  over  the  portion  of 
the  time  history  in  which  q(t)  increases  monotonically .  However,  it  will  be 
found  that  this  interval  encompasses  not  only  the  interval  over  which  tum¬ 
bling  occurs,  but  also  the  subsequent  range  over  which  the  oscillatory  motion 
begins  and  is  reduced  to  small  angles.  The  results  of  this  analysis  should 
be  suited  to  act  as  the  connecting  link  between  the  bodyfs  initial  behavior 
and  the  behavior  described  by  the  results  of  analyses  in  which  the  linear 
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approximation  to  the  aerodynamic  restoring  moment  has  been  introduced  (ref.  5; 
e.g.,  which  is  applicable  after  the  oscillatory  motion  has  been  reduced  to 
small  angles). 


Transformed  equation. 

gives 


Inserting  equations  (5)  and  (8)  into  equation  (4) 


Al  n  Jv±t 
^i  i  ^“max*5 


sin  a  =  0 


which  will  be  taken  to  be  the  differential  equation  characterizing  the  turn- 
bling  entry  problem.  However,  a  transformation  of  equation  (9)  yields  a  form 
that  more  quickly  shows  the  nature  of  the  solution.  Let 


PVi  =  s 

Mr  -  ( lliY 

-(^i  I  “max  \  2 )  ) 


X 

k 


eSt/2  J 


(10) 


Equation  (9) 


takes  the  form 


a"(x)  + 


a*(x)  + 

X 


sin  a  =  0 


(ii) 


with  the  initial  conditions 


a(  k  )  = 

a'(K)  =  h 


(12) 


where,  for  convenience,  a(ic)  is  presumed  to  lie  within  the  range  ~k  <  a(ic)  <  rt. 
Note  that  all  the  parameters  of  the  problem  have  been  concentrated  in  the  con- 
stant  K  and  the  initial  conditions.  Equations  (ll)  and  (12)  indicate  that 
all  combinations  of  body  and  planetary  properties  yielding  the  same  value  of 
k  and  the  same  initial  conditions  og(k),  a' (k)  will  yield  identical  solutions 
for  a  as  a  function  of  x,  though  not  necessarily  as  a  function  of  time. 

The  independent  variable  x  may  be  given  a  physical  meaning  when  it  is  noted 
that,  within  the  approximations  made,  it  is  proportional  to  the  square  root  of 
the  dynamic  pressure  q.  .  . 


The  Painleve  Transcendents 


The  substitution  w  =  sin  a/2  in  equation  (ll)  transforms  it  to 

w"(x)  =  L(w)p2  +  M(x)p  +  N(w)  (13) 
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where 


L(w) 


w 


W2  -  1 


M(x)  =  -  i 
N(w)  =  w(w2  -  1) 


In  reference  6^  it  will  be  found  that  equation  (13)  is  of  the  form  studied  by 
the  French  mathematician  Paul  Painleve  around  1900.  It  is  categorized  mathe¬ 
matical]^-  as  being  a  member  of  the  general  class  of  second-order  differential 
equations  whose  solutions  have  fixed  critical  points  (i.e.,  no  movable  branch 
points  or  essential  singularities).  This  class  has  been  shown  to  possess  50 
members.  Of  the  50,  all  but  6  are  integrable  in  terms  of  known  functions. 

The  remaining  6  define  new  functions,  termed  "Painleve  transcendents."  The 
substitution  W  =  (w  +  l)/(w  -  1)  casts  equation  (13)  in  the  form 


W"(x) 


w  +  i\ 

W  -  l) 


(1^) 


1  rng 

Sketch  (b) 


and  it  will  be  seen  in  reference  6 
that  equation  (l4)  is  one  of  these, 
namely,  the  fifth  (with,  in  Ince!s 
notation,  a  =  3  =  7  =  0,  b  =  -2). 
Unfortunately,  aside  from  this  categor¬ 
ization  and  a  comprehensive  examination 
of  the  asymptotic  behavior  of  the  first 
Painleve  transcendent  (ref.  7 ) >  n0  sub¬ 
sequent  analyses  of  their  properties 
seem  to  have  been  published. 


Mechanical  Analogy 

Before  proceeding,  let  us  first 
consider  a  simple  mechanical  analogy  of 
equation  (9)  whose  behavior  is,  in 
effect,  intuitively  obvious.  By  this 
means,  the  range  and  character  of 
motion  governed  by  the  equation  can  be 
revealed  relatively  simply. 

Consider  a  small  bead  constrained 
to  slide,  without  friction,  on  a  circu¬ 
lar  path  in  a  vertical  plane.  Let  a 
time -dependent  force  F(t)  be  exerted 
downward  on  the  bead.  This  situation 
is  illustrated  in  sketch  (b).  Equating 
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■the  torque  about  the  center  of  the  circle  to  the  bead’s  rate  of  change  of 
angular  momentum  gives 
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1 

+  — 
r 


g  + 


F(t) 


m 


sin  cp  =  0 


(15) 


Thus,  if  g  +  [F(t)/m]  is  caused  to  vary  in  proportion  to  the  body's  dynamic 
pressure  history,  the  motion  cp(t)  governed  by  equation  (15)  will  be  analogous 
to  that  of  the  body  with  9  playing  the  role  of  a. 

For  convenience,  let  cp  =  0  (the  bottom  of  the  circle)  at  t  =  0.  Now 
give  the  bead  an  initial  velocity  sufficient  to  carry  it  several  times  around 
the  circle.  Each  time  the  bead  traverses  the  circle,  cp  will  be  counted  as 
having  increased  by  2rt .  This,  of  course,  corresponds  to  one  complete  tumble 
of  the  body.  It  is  clear  that,  because  F(t)  increases  continually,  more  and 
more  of  the  total  energy  will  be  in  the  form  of  potential  energy  each  time  the 
bead  nears  the  top  of  the  circle.  Eventually,  therefore,  the  bead  will  not 
have  sufficient  kinetic  energy  to  carry  it  over  the  top.  At  this  point,  it 
will  reverse  direction,  slide  down  past  the  low  point,  and  proceed  to  oscil¬ 
late  about  the  low  point.  Again,  because  the  amplitude  of  the  restoring, 
torque  grows  indefinitely  with  time,  the  amplitude  of  oscillation  will  dimin¬ 
ish  and  the  frequency  will  increase  with  time.  The  final  value  of  cp  will  be 
a  multiple  of  2jt .  This  behavior  is  illustrated  as  curve  A  in  sketch  (c)  for 
a  case  where  the  bead  has  tumbled  twice. 

Note  in  the  sketch  that  once  tumbling  is 
arrested,  the  bead's  amplitude  of  oscil¬ 
lation  about  2nrt  cannot  exceed  jc.1 
For  a  range  of  successively  smaller  ini¬ 
tial  velocities,  the  behavior  of  cp(t) 
will  be  qualitatively  similar  to  that 
just  described,  the  tumbling  in  each 
case  being  arrested  when  cp  is  between 
(2n  -  l)jt  and  (2n  +  l)jt  and  the  subse¬ 
quent  oscillation  being  about  2mt . 

Eventually,  however,  as  the  initial 
velocity  is  successively  reduced,  a 
specific  initial  velocity  will  be 
reached  for  which  the  kinetic  energy 

near  the  top  of  the  circle  (i.e.,  cp  =  (2n  -  l)jt)  is  just  sufficient  to. enable 
the  bead  to  reach  the  top  and  come  to  rest  there.  As  the  top  is  a  position 
of  unstable  equilibrium,  the  bead  cannot  oscillate  about  that  position,  but 
must  approach  it  uniformly  from  below.  This  is  shown  as  curve  B  on 
sketch  (c).  For  a  slightly  smaller  initial  velocity,  the  bead  will  not _ 

ijt  will  be  observed  that  this  description  of  the  arrest  of  tumbling  does 
not  entail  the  presence  of  aerodynamic  damping.  This  contradicts  a  result 
presented  in  reference  8,  in  which  the  cause  of  the  arrest  of  tumbling  is 
attributed  to  a  dissipation  of  rotational  energy  through  aerodynamic  damping. 
The  author  of  reference  8  is  led  to  this  conclusion  by  the  erroneous  assump¬ 
tion  that  the  net  change  of  potential  energy  over  one  complete  revolution  is 
small  enough  to  be  neglected. 
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Sketch  (c) 
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surmount  the  top  and,  hence ,  must  oscillate  about  the  next  smaller  multiple 
of  2rt  (i.e.,  (2n  -  2)jt).  This  is  shown  as  curve  C  on  sketch  (c)  .  Again, 
there  will  be  a  range  of  successively  smaller  initial  velocities  for  which 
the  bead  will  oscillate  about  (2n  -  2)jt,  eventually  terminating  with  an  ini¬ 
tial  velocity  that  brings  the  bead  to  rest,  without  oscillating,  at  (2n  -  3)* 
(curve  D) . 

To  summarize  this  discussion  as  it  applies  now  to  the  body,  one  observes 
that  for  any  initial  angle  of  attack  there  will  be  a  range  of  values  of 

initial  angular  velocity  dq  that  will  cause  oo  eventually  to  oscillate 
about  a  given  value  of  2mt .  This  range  of  initial  angular  velocities  is 
bounded  by  the  two  specific  angular  velocities  which,  for  the  same  aq,  cause 
a  to  come  to  rest  without  oscillating  at  (2n  +  l)rt  and  (2n  -  l)jt.  The  latter 
cases,  where  the  body  comes  to  rest  in  a  position  of  unstable  equilibrium,  are 
somewhat  unrealistic  in  practice.  They  can  be  important  in  analysis,  however, 
as  they  serve  to  define  the  multiple  of  2-x  about  which  a  eventually 
oscillates. 


Properties  of  the  Transcendent 


Although  it  is  possible  to  obtain  approximate  analytical  representations 
of  the  motions  just  described  (cf.  ref.  2),  these  solutions  cannot  be  extended 
far  enough  to  cover  all  cases  of  interest.  In  particular,  they  begin  to  fail 
as  they  approach  the  condition  for  which  the  body  tends  to  dwell  near  a  posi¬ 
tion  of  unstable  equilibrium.  On  the  other  hand,  it  is  not  useful  merely  to 
present  a  catalog  of  numerical  solutions  for  the  Painleve  transcendent,  as  a 
separate  solution  would  be  necessary  for  each  pair  of  initial  conditions  and 
this  would  require  a  prohibitive  number  of  solutions.  Alternatively,  analysis 
can  point  the  way  to  obtain  a  limited  number  of  generalized  numerical  results 
from  which  most  of  the  properties  of  interest  can  be  derived.  This  will  be  ' 
the  objective  of  the  present  section. 

Definition  of  parameters.  -  Consider  equation  (ll),  and  to  bring  in 
evidence  that  for  the  oscillatory  solution,  a  ->  2nrt  as  x  -»  00 ,  let 

oc  =  2njt  -  2e  (l 6) 

where  n  is  an  integer,  denoting  the  number  of  tumbles  completed  over  the 
range  of  tumbling.  As  indicated  in  the  preceding  discussion,  if  its  value  is 
required,  it  must  be  determined  from  the  nonos dilatory  solution  (cf.  ref.  2). 
In  the  present  study,  however,  it  will  not  be  necessary  to  know  the  cumulative 
value  of  a,  but  only  its  value  relative  to  the  -final  equilibrium  position 
2m*  Hence,  n  need  not  be  determined.  It  is  convenient  to  retain  the  symbol 
a  for  designating  angle  of  attack,  but  in  doing  so,  henceforth  we  shall 
actually  mean  a  -  2n^:,  or  interchangeably,  -2e.  The  equation  of  motion  in 
terms  of  e  becomes 

e”  +  —  +  sin  e  cos  e  =  0  (17) 

x 
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Locus  of  first  peaks,  x=x(£p) 


mm 


where  now  e  ->  0  as  x  -*■  Consider  a 
sequence  of  solutions  of  equation  (17)  , 

in  which  the  initial  velocity  €*(k)  is 
held  constant  and  e(rc)  is  allowed  to  -7r/2 
vary.  A  typical  sequence  is  shown  on 
sketch  (d).  Observe  that  for  the  given 
€'(k),  there  is  a  unique  value  of  e(ft) 
within  an -interval  of  it  that  will 
bring  the  motion  to  rest  without  oscil-  T/2 
lating  at  e  =  -it/2  (corresponding  in 
a  to  an  odd  multiple  of  it).  This  €  e 

motion  is  shown  as  curve  A  in  sketch  (d)^ 
and  its  initial  value  will  be  called  the 
critical  angle,,  ec(k).  It  should  be 
clear  that  the  curve  labelled  B  is  iden¬ 
tical  to  curve  A  and  describes  the  same 
event,  as  do  any  two  curves  in  e  whose  «c+7r 
initial  velocities  €'(k)  are  the  same  x=k 

and  whose  initial  values  are  separated 

by  a  multiple  of  it.  As  is  evident  from  Sketch  (a) 

the  sketch,  all  possible  solutions  for  a 
given  e 1 ( k )  are  included  within  and 

bounded  by  the  nonoscillatory  solutions  shown  as  curves  A  and  B.  Now  in- most 
applications,  the  questions  of  greatest  interest  will  be:  When  does  tumbling 
effectively  stop  and  what  is  the  nature  of  the  subsequent  oscillation?  Both 
questions  can  be  answered  by  defining  conditions  at  the  first  peak  of  the 
oscillatory,  motion.2  As  shown  on  sketch  (d),  a  locus  of  such  points  can  be 
drawn  which  then  intersects  the  first  and  largest  peak  of  every  one  of  the 
oscillatory  motions.  For  any  s(k)  between  ec  and  ec  +  it ,  let  €p  be  the 
value  of  e.  at  the  first  peak  and  x(ep)  be  the  value-  of  x  at  the  first 
peak.  Functionally,  ep  and  x(ep)  are  dependent  on  the  initial  conditions; 

th.Eit  IS  .  f  f  \  1  /  ,  \  \  1  • 

€p  =  ep(e(ft)>Ke  (k)jk)  7 

■  *  p  (18) 

x(ep)  =  x(e(K),Ke'(K),K) 

For  entry  from  without  the  atmosphere,  for  which  k  -*■  0,  it  can  be  shown  that 
the  functional  relationships  (l8)  no  longer  depend  on  K.  Further,  it  will  be 
found  convenient  to  show  the  dependence  On  e(ic)  as  a  dependence  on  an  incre¬ 
ment  in.  e ( k ) ,  the  difference  between  e(K)  and  the  critical  angle  ec.  Thus, 

let  ■  : 

Ae.i  =  e(fe)  -  ec  • ,  .  (19) 

For  entry  from  without  the  atmosphere,  then 

.  ,  =  ep(u-i/s,Aei)  I  ,  (20) 

■  x(ep)  =  x(aj/s,Aei)  j 

be  precise,  tumbling  actually  stops  earlier  than  at  the  first  peak,  - 
namely,  at  the  value  of  x  for  which  e,2(x)  =  cos2  e  (cf.  ref.  2).  However, 
the  difference  between  the  two  coordinates  generally  is  not  great  and  the 
location  of  the  first  peak  is  a  physically  more  tangible  parameter  to  charac- 
terize  the  end  of  tumbling. 
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since 


Ke1  (k)  =  -cy_/s 

Note  that  these  functional  relationships  are  independent  of  all  body  param¬ 
eters.  They  indicate  the  existence  of  single  solutions  applicable  to  all 
cases.  Thus,  it  is  practical  to  evaluate  them  numerically  once  and  for  all. 
This  has  been  done  and  the  results  are  presented  on  figures  1  and  2.  In  an 
entry  from  without  the  atmosphere,  all  values  of  Ae±  are  equally  probable; 
however,  the  values  of  primary  concern  are  those  which  lead  to  the  largest 
first  peaks  of  oscillation,  for  the  subsequent  oscillation  amplitudes  will 


-90 


-50 


-40 


Figure  1.-  Values  of  e  at  the  first  peak  of  oscillation  as  a  function  of  initial  angular 

velocity. 
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28 


Figure  2.-  Values  of  x  at  the  first  peak  of  oscillation  as  a  function  of  initial  angular 

velocity. 


then  be  largest  through  the  high  heating  and  loading  portions  of  the  trajec¬ 
tory.  Reference  to  sketch  (d)  will  show  that  the  largest  first  peaks  occur 
when  Ae±  is  near  zero  and  again  near  «,  and  values  of  ep  and  x(ep)  for 
these  ranges  of  Aej_  have  been  particularly  stressed  in  figures  1  and  2.  It 
will  be  seen  also  on  sketch  (d)  that  there  exist  two  solutions  of  special  sig¬ 
nificance.  One  of  these  gives  the  minimum  possible  value  of  x  for  which  the 
first  peak  can  occur.  It  is  denoted  by  xPfflin  and  e(xpmin)  on  figures  1  and 
2.  The  second  special  solution  leads  to  the  smallest  possible  value  of  e  at 
the  first  peak.  It  is  denoted  by  epmin  and  x(ePmin)  on  1  and  2>  The 

latter  solution  is  particularly  important  since,  as  will  be  shown,  it  may  be 
used  to  establish  a  lower  bound  for  the  envelope  of  oscillatory  motion.  That 
this  minimum  angle  of  attack  at  the  first  peak  of  oscillatory  motion  must  be 
different  from  zero  will  be  evident  when  it  is  recalled  that  all  of  the  rota¬ 
tional  energy  of  the  body  must  show  up  as  potential  energy  due  to  angle  of 
attack  when  the  body  pitch  rate  is  zero. 

Asymptotic  solution  for  zero  initial  pitch  rate. -  With  the  magnitude  and 
coordinate  of  the  first  peak  of  oscillatory  motion  defined,  it  is  now  possible 
to  treat  the  oscillatory  motion  relatively  simply.  In  effect,  the  problem  has 
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been  reduced  to  the  case  of  a  nontumbling  body  which  begins  its  oscillatory 
motion  at  x(ep),  a  particular  value  of  x,  with  the  initial  conditions 
€i  =  6p,  e±*  =  0.  Let  f  =  tan  e.  Equation  (17)  becomes 


p-p-p  r  2  .pr 

fTT-  — — 5-  +  —  +  f  =  0  (21 

1  +  f 2  x 

Noting  the  resemblance  of  equation  (21)  to  the  equation  for  the  Bessel  func¬ 
tion  of  zero  order,  one  may  cast  equation  (21)  in  the  form  of  an  integral 
equation  involving  the  Bessel  function  by  means  of  the  method  of  variation  of 
parameters .  Thus , 

pOO 

f  =  Z0(x)  +  |  J  tl( 5)[J0(x)Y0(i)  -  Y0(x)Jo(|)]d|  (22; 


with 


z(x)  =  I 

1  +  f2(x) 

Z0(x)  =  U0(x)  +  bY0(x) 


pOO 

-|  J  |2(|)Y0(|)d| 

x(£p) 


b  =  b  +  | 


x(ep) 


a  =  -  tt  x 


( £-d)Yt  (x( e-n) )  tan 


2  \  ^p  /  x x  n  ^p 


b  =  -  x(ep)J1(x(ep) )  tan  ep 


The  integral  in  equation  (22)  can  be  shown  to  diminish  faster  than  x~1/2  as 
x  ->  °o,  and  hence  the  asymptotic  behavior  of  f  (x)  is  f (x)  ~  Z0(x) .  Values 
of  x  in  the  portion  of  the  trajectory  of  primary  interest,  where  aerody¬ 
namic  loads  and  heating  become  significant,  are  sufficiently  large  that  the 
asymptotic  behavior  of  ZQ(x)  may  be  used.  Then,  for  large  x 


^  (a2  +  b2)  sin  -  j-  + 


where 


tan'-1  E1 
b 


The  envelope  of  oscillation  is  given  by 


(tan  e)gnv  ~  (a  +  ^  ) 


(25) 


Since  for  small  values  of 
equation  (25)  as  follows: 


:P' 


a  ->  a 


b  b,  it  is  convenient  to  normalize 


(tan  e)env 


J 


; (a2  +  b2) 


a2  +  b2 
a2  +  b2 


(26) 


Let 


where 


and 


where 


G 


ag  +  b5 
a2  +  b2 


G  =  G(ep,x(ep)  ) 


(27) 


c  =  x(ep) J |  [Yi2(x(ep))  +  Ji2(x(ep))] 


C  =  C(x(ep)) 

Then  equation  (26)  takes  the  simple  form 

4x  (tan  e) 


J 


env 


=  G 


C  tan  e- 


(28) 


Equation  (28)  is  the  correct  asymptotic  expression  for  the  envelope  curve  of 
tan  e  in  terms  of  the  parameter  G.  Unfort unately,  G  cannot  be  evaluated 
analytically  but  since  it  is  dependent  upon  only  two  parameters,  ep  and  x(ep), 
it  is  practical  to  evaluate  it  numerically  for  a  wide  enough  range  to  suffice 
for  all  conditions  of  interest.  This  has  been  done  and  the  results  axe  pre¬ 
sented  in  figure  3. 3  For  convenience,  the  parameter  C  has  also  been  evalu- 
ated  and  the  results  are  presented  in  figure  t.  With  G  determined,  equation 
(28)  is  a  very  general  expression  for  the  body's  envelope  of  oscillatory 
motion.  Observe  that  it  is  independent  of  all  body  parameters.  It  is  spe¬ 
cialized  to  a  particular  case  by  relating  x  to  some  physically  significant 
parameter  such  as  altitude,  fraction  of  maximum  dynamic  pressure,  or  fraction 
of  maximum  aerodynamic  heating.  Examples  of  these  relationships  will  be  given 
in  the  next  section. 


3The  origins  of  the  curves  for  x(ep)  =  10  and  x(ep)  -  20  have  been 
displaced  for  clarity. 
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Finally,  consider  the  case  x(ep)  ->  0,  which  corresponds  to  the 
nontumbling  entry  of  a  body  from  without  the  atmosphere.  Inspection  of  the 
results  of  figure  3  shows  G  to  be  independent  of  x(ep)  for  values  of  x(ep) 
approaching  zero.  Further,  C  (eqs.  (27))  is  also  independent  of  x(ep)  for 
small  x(ep),  as  may  be  shown  when  the  Bessel  functions  in  C  are  replaced 
with  their  initial  behavior  approximations.  Thus, 

C  -  for  x(€p)  -  0  (29) 

With  these  results,  equation  (28)  may  be  further  simplified  for  the  case  of  a 
nontumbling  entry  from  without  the  atmosphere.  The  envelope  of  oscillatory 
motion  is  given  by 
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(30) 


(tan  e)env  =  G(ep)  tan  eP 


It  is  of  particular  interest  to  note  the 
general  nature  of  the  solution  offered 
in  equation  (30).  The  solution  written 
as  a  function  of  the  coordinate  x  is 
independent  of  all  body  parameters, 
planetary  properties,  and  initial  condi¬ 
tions  except  initial  angle  of  attack. 


Useful  Relations 


Within  the  assumptions  of  this 
analysis  (V  =  V±,7  =  7±) >  the  dynamic 
pressure  varies  as 

q/q.  =  ePVi  Sln  ^  (3l) 

Hence,  from  equations  (10) 


0  4  8  12  16  20 


X<€p) 

Figure  4.-  Variation  of  the  parameter  C 
with  x(ep) . 


X 

K 


so  that  the  envelope  expression  (eq.  (28))  may  be  rewritten 


(32) 


(tan  Oenv  =  CGk-i/2  M 
tan  €p  \qi/ 


(33) 


It  will  be  noted  that  this  form  is  in  agreement  with  the  asymptotic  form  (for 
negligible  aerodynamic  damping)  previously  derived  in  reference  5  for  the 
case  of  small  initial  a.  The  latter  derivation  is  not  dependent  on  assump¬ 
tions  concerning  the  variations  of  either  flight -path  angle  or  flight  veloc¬ 
ity.  Hence,  equation  (33)  can  be  extended,  and  the  restrictions  underlying 
its  development  circumvented,  by  the  simple  expedient  of  using  a  more  precise 
expression  for  q/q.  in  place  of  equation  (31)-  The  results  of  reference  9 
are  particularly  useful  in  this  regard  as  they  permit  one  to  write  this  ratio 
in  a  variety  of  forms  involving  physically  significant  parameters. 


Envelope  of  oscillation  in  terms  of  altitude. -  Erom  reference  9 


\ 

so  that 


V  =  V^e 
2 


-(CpP0A/2(Bm  sin  7j:)e 


•Py 


0~V-r  ( n^r\  A  /Rtvi  -?  m 


-3y 


(34) 


Substituting  equation  (35)  in  (33)  gives 


where 


(tan  e)eny 
tan  e-n 


=  CG 


sin2  7il  rPyp  (CpP^A/Pm  sin  y.  )e 
L  2p0AZC; 


-Py 


-,1/4 


Lmax 


(36) 


€P  =  €p(°c,i/s^^€i)  y  C  -  C(x(ep))  ,  G  -  G(ep;x(€p)) 


It  should  be  kept  in  mind  in  using  equation  (36)  that  x  has  been  assumed  to 
be  sufficiently  large  that  the  solution  is  indeed  described  by  its  asymptotic 
form.  Thus,  there  is  implied  an  upper  bound  on  y  which  will  exclude  values 
of  y  large  enough  to  cause  (tan  e)env/tan  Gp  to  exceed  or  approach  unity. 


Envelope  of  oscillation  in  terms  of  dynamic  pressure. -  From  reference  9 
the  maximum  dynamic  pressure  is 


Writing  q/q^  as 


i  (Vax 


1 

2e 


-i— -) 

sin  7. -J 


JL  = 

*1 


q  ^max 
%nax 


(37) 

(38) 


and  substituting  equations  (37)  and  (38)  in  (33)  gives 


(tan  e)eny. 
tan 


=  CG 


1 

2  VA2Cmi 


(3e  sin  y±) 


'CpAN  a 


m 


-,1/4 


Lmax 

q. 


(39) 


Envelope  of  oscillation  in  terms  of  aerodynamic  heating*-  The  results  of 
reference  9  provide  the  following  expression  for  fraction  of  maximum 
stagnation -point  heating  rate 


js 

^s^max 


-Ptor/2)  (i/2)(!-e-p  Ay) 

e  “ 


where 

=  y  -  yi 

Yx  =  altitude  for  maximum  heating  rate  S 
_  1 

yi  P  \pm  sin  jtJ  - 


(^0) 


(41) 
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Equation  (1+0)  cannot  be  inverted  to  solve  for  £y,  so  that  a  graphical  solu¬ 
tion  is  necessary.  This  is  presented  in  reference  9.  Substituting  equa¬ 
tions  (1+1)  in  (36)  gives 


(tan  e)gnv 
tan  ep 


CG 


.  i  sin  / _ I _ V^e(i/3)e‘P^ll/4 

2Ur  s'n7iVCBta/  J 


where  0  Ay  can  be  related  to  Hs/(Hs)max  through  the  graphical  solution 
provided  in  reference  9« 


APPLICATION’S 


As  illustration  of  their  use,  the  preceding  results  will  be  applied  in 
two  different  connections;  first,  as  an  aid  in  the  design  of  vehicles  intended 
for  use  as  planetary  probes  and,  second,  as  an  adjunct  to  the  research  into 
the  origin  of  tektites.  As  a  preliminary  to  both  applications,  however,  a 
method  will  be  presented  that  will  be  of  general  use  in  placing  reasonable 
bounds  on  the  oscillation  envelopes  of  initially  tumbling  bodies. 


Bounds  on  Oscillation  Envelopes 


A  body  making  an  uncontrolled  entry  into  a  planetary  atmosphere 
generally  will  be  tumbling  at  an  essentially  constant  rate  prior  to  entering 
the  sensible  atmosphere.  Because  of  the  tumbling  motion,  generally  it  will 
not  be  possible  to  specify  the  angle  of  attack  at  the  precise  moment  aerody¬ 
namic  effects  begin  to  influence  the  motion.  Therefore,  an  uncertainty  in  the 
angle -of -attack  history  within  the  atmosphere  will  necessarily  prevail.  How¬ 
ever,  it  is  possible  to  relate  the  parameters  describing  the  body  and  its 
trajectory  to  the  probability  that  the  body  will  eventually  undergo  oscilla¬ 
tions  bounded  by  specified  envelopes.  Thus,  the  uncertainty  in  the  angle-of- 
attack  envelope  history  within  a  given  atmosphere  is  bounded  below  by  a 
minimum  possible  for  a  given  initial  tumbling  rate  and  above  by  a  maximum  for 
which  the  probability  of  exceeding  the  maximum  can  be  specified. 

Consider  first  the  lower  bound.  The  results  giving  envelopes  of 
oscillation  in  terms  of  aerodynamic  heating  and  dynamic  pressure  may  be  put 
into  a  particularly  convenient  form.  First,  let 


B 


(|B  sin 


7  i) 


AZCm 


■max 


m 


Then,  with  the  use  of  equations  (39)  and  (1+2),  the  values  of  the  minimum  pos¬ 
sible  envelopes  evaluated  at,  respectively,  maximum  dynamic  pressure  and 
maximum  heating  may  be  given  the  form 
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(tan  e)m^n  env 
(tan  e). 


max  q 


rain  env 


i/4 

=  GG[  -  —  B  )  tan  e 
.  1/4 

_  3  .1/3 


ram 


=  CG  -  —  e'L/  B  tan  e 


•max  heat 


ram 


m 


Note  that,  aside  from  the  parameter  B,  equations  (44)  contain  terms  involving 
only  the  initial  tumbling -rate  parameter  2dj_/s.  Thus,  if  values  of  B  are 
chosen  for  a  wide  enough  range  to  cover  all  cases  of  interest,  these  relations 
may  be  evaluated  once  and  for  all.  The  results  are  presented  on  figure  5* 

Now,  it  is  further  noted  that  the  following  relations  hold 


(tan  e) 


mm  env 


(tan  e)  • 


min  env 


Lmax 

q 


1/4 


max  q 


(tan  o)in^n  env 


(tan  e) 


rain  env 


max  heat 


P  Ay  (i/3)e“^  ^ 


e  v  e 


1/3 


1/4 


>  (45) 


and  that  these  expressions  are  independent  of  all  parameters  relating  to  the 
body.  These  functions  are  presented  graphically  in  figure  6.  The  minimum 
possible  envelope  of  oscillation  in  terms  of  fraction  of  maximum  dynamic 
pressure  or  fraction  of  maximum  heating  can  be  computed  from  the  graphs  of 
figures  5  and  6.  A  relationship  between  the  minimum  envelopes  and  envelopes 
having  a  specified  probability  of  being  exceeded  will  now  be  established. 


A  study  of  sketch  (d),  for  a  given  initial  tumbling  rate,  suggests  the 
possibility  that  the  asymptotic  behavior  of  the  envelopes  of  oscillatory 
motion  for  two  different  initial  angles  of  attack  might  be  identical. 
Inspection  of  equation  (28)  reveals  that  this  will  indeed  be  the  case  if  the 
term  CG  tan  ep  is  the  same  for  two  different  initial  angles  of  attack.  Now 

it  is  noted  that  this  product  is  a  function  only  of  Ae±  for  a  given  tum¬ 
bling  rate  and  that  for  each  value  it  takes  on  for  Aei  near  zero  there  is  a 

corresponding  Aei  near  Jt  for  which  it  has  the  same  value.  Values  of  Aei 

near  zero  correspond  physically  to  the  body  being  oriented  in  an  attitude 
slightly  less  inclined  to  the  stream  than  the  critical  attitude  and  values  of 
Aei  near  tc  correspond  to  the  body  being  at  a  slightly  greater  inclination 
than  for  the  critical;  an  appropriate  interpretation  of  these  values  gives 
the  probability  of  exceeding  a  particular  envelope.  Thus,  suppose  that  a 
particular  value  of  CG  tan  ep  is  chosen  and  that  this  same  value  results 
from  both  Aei  =  1*5°  and  Aei  =  179 • 5° •  For  values  of  Aei  between  0°  and 
1.5°,  or  between  I79.50  and  l80°,  the  subsequent  envelopes  of  oscillation 
will  exceed  those  for  all  Aei  between  1*5°  and  179* 5°*  Then  there  is  a 
2.0°  range  of  Aei  °ut  of  a  possible  l80°  range,  or  a  probability  of  l/90 
that  the  envelope  for  the  chosen  value  of  CG  tan  ep  will  be  exceeded. 
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Figure  6.-  Relationship  between  minimum  possible  envelope  values  of  tan  e  and  fraction 
of  maximum  heating  and  maximum  dynamic  pressure. 


The  relationship  between  an  envelope  of  oscillation  having  a  specified 
probability  of  being  exceeded  and  the  minimum  possible  envelope  is  independ¬ 
ent  of  the  characteristics  of  the  body.  One  observes  this  by  forming  the 
following  ratio,  using  equation  (28) 

C(ep)G( €p,x( Sp) )  tan  Gp 

CGPmin}GC€Pmin?^ePmin))  tan  €pmin 

and  noting  that  the  right  side  is  dependent  only  upon  the  tumbling  rate 
parameter  2dp/s.  This  relationship  has  been  evaluated  for  several  probabili¬ 
ties  and  the  results  are  presented  on  figure  7*  Application  of  the  results  of 
figures  5*  6,  and  7  to  make  estimates  of  the  upper  and  lower  bounds  to  the 
envelope  of  oscillation  is  illustrated  in  the  following  section. 


(tan  e)env 
an  G)min  env 
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(tan  Oenv  /  (tan  €)min  env 


2dj/s 


Figure  7.-  Relationship  between  envelope  values  of  tan  e  having  a  specified  probability  of 
being  exceeded  and  values  of  tan  e  for  the  minimum  possible  envelope  as  a  function  of 
initial  angular  velocity. 


Motions  of  Planetary  Probe  Vehicles 


To  illustrate  the  use  of  the  results  ana  to  assess  their  accuracy,  a 
sample  calculation  will  be  made  ana  compared  with  results  in  which  no  approxi¬ 
mations  were  made  either  to  the  equations  of  motion  or  to  the  aerodynamic 
forces  and  moments.  The  comparison  will  be  made  with  the  results  of  refer¬ 
ence  1  for  a  vehicle  making  a  tumbling  entry  into  the  martian  atmosphere. 

Calculation  of  parameters. -  The  vehicle  studied  in  reference  1  had  the 
following  physical  properties: 


A  =  8.296  ft2 

I  =  5*6  lb -ft -sec2 


l  =  3.25  ft 
m  =  6.685  lb-sec2/ft 

cp  =  0.650 


C 


'■''max 


-O.162 


where  Cmrnay  was  evaluated  so  that  the  area  under  a  half -cycle  of  the  approx¬ 
imating  sine  function  curve  equalled  that  under  the  actual  pitching -moment 
variation.  The  martian  atmosphere  and  the  entry  conditions  were  given  as 
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leading  to 


0  =  2 . 15x10  ~5/ ft 
V±  =  21,042  ft/sec 

7i  =  41-5° 

=  0.2094  radian/sec 


s  =  0.30/sec 
2di/s  =  1.396 

B  =  -1.43 2x10 "5 
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Figure  8.-  Envelopes  of  oscillatory  motion 
as  a  function  of  fraction  of  maximum  heat¬ 
ing  rate  for  a  Mars  entry;  initial  angular 
velocity  =  12°/sec. 


Determination  of  angle -of -attack 
envelopes . -  The  minimum  possible  enve¬ 
lope  of  oscillatory  motion  in  terms  of 
aerodynamic  heating  rate,  subsequent  to 
the  arrest  of  tumbling,  is  jdeter  mined 
in  the  following  way:  Enter  figure  5 
with  the  given  values  of  2ai/s  and  B 
to  obtain  [(tan  e)min  envlmax  heat* 

Use  this  value  in  conjunction  with  fig¬ 
ure  6  to  construct  a  graph  of  (a)min  env 
versus  Hs/HSmax.  Results  of  this  cal¬ 
culation  are  presented  on  figure  8  and 
the  exact  results  are  shown  for  compar¬ 
ison.  One  obtains  the  envelope  having 
a  probability  of  2.75/360  of  being 
exceeded  by  entering  figure  7  with  the 
given  values  of  2di/s  and  probability; 
this  gives  the  ratio  of  the  envelope 


values  of  the  upper  and  lower  bounds  of  tan  e.  The  upper  bound  envelope  is 
constructed  from  the  product  of  this  ratio  and  the  already  obtained  minimum 
envelope.  Results  of  this  calculation  are  also  shown  on  figure  8  and  compared 
with  the  exact  results .  It  is  noted  that  the  agreement  between  the  exact 
result  and  that  of  the  present  theory  is  within  10  percent  throughout. 


Motions  of  Tektites 


Origin  of  tektites. -  From  the  exhaustive  studies  of  the  nature  of  the 
curious  glass  objects  known  as  tektites  (see  ref.  4  for  an  expensive  bibliog¬ 
raphy)  ,  it  has  been  known  for  a  long  time  that  the  tektites  experienced  two 
separate  periods  of  intense  heating.  During  the  first,  heating  was  suffi¬ 
ciently  intense  to  melt  the  objects  completely,  whereas  in  the  second,  heating 
was  sufficient  only  to  melt  the  thin  surface  layers  of  otherwise  solid  objects. 
After  the  extensive  investigations  of  Chapman  and  his  colleagues  (refs.  4  and 
10),  in  which  the  ring-wave  markings  and  coiled  flanges  that  are  such  distinc¬ 
tive  features  of  the  Australian  tektites  were  reproduced  with  fidelity  in 
wind-tunnel  experiments,  there  can  be  little  doubt  that  the  second  period  of 
melting  occurred  as  the  result  of  a  hypervelocity  passage  through  the  Earth* s 
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atmosphere.  Further,  from  a  minute  study  of  the  striae  lying  beneath  the 
tektites ’  surfaces,  combined  with  the  observed  final  tektite  geometry,  the 
authors  of  references  4  and  10  have  been  able  to  deduce  the  probable  speed  and 
flight -path  angle  of  these  objects  as  they  entered  the  atmosphere,  and  they 
have  concluded  that  the  entry  conditions  deduced  were  compatible  with  those 
for  objects  whose  origin  was  the  Moon.  If  this  conclusion  is  accepted. as  cor¬ 
rect,  the  following  sequence  of  events  may  be  hypothesized:  The  Moon. is 
struck  by  a  meteor,  the  lunar  surface  at  the  point  of  impact  is  vaporized  and 
streams  of  molten  lunar  surface  material  are  ejected  outward  from  the  crater 
at  very  great  velocity,  at  such  great  velocity  indeed  that  the  molten  material 
has  sufficient  kinetic  energy  to  escape  the  gravitational  field  of  the  Moon. 
Subsequently,  some  of  the  lunar  material  traverses  a  path  in  space  that  brings 
it  within  a  corridor  permitting  capture  by  the  Earth’s  gravitational  field. 
Sometime  after  its  ejection  from  the  Moon,  the  stream  of  molten  material 
breaks  into  segments  which  tend  to  contract  into  characteristic  shapes  by  the 
action  of  surface  tension,  congeal,  and  solidify  as  they  lose  their  heat  by 
radiation.  The  tektites  then  enter  the  Earth's  atmosphere  as  solid  bodies, 
and,  during  their  passage  through  the  atmosphere,  acquire  by  ablative  melting 
the  characteristic  ring  waves  and  flanges  that  so  distinguish  them  when  they 
are  found  on  the  Earth’s  surface. 

Now,  if  this  sequence  did  actually  originate  with  a  collision,  it  seems 
probable  that  at  least  a  portion  of  the  lunar  material  should  have  acquired  a 
certain  amount  of  angular  momentum.  The  molten  objects  then  would  have  found 
themselves  turning  while  in  space  and  would  have  sought  to  assume  the  stable 
figures  consistent  with  their  turning  rates.  They  would  have  been  turning, 
or  tumbling,  also  as  they  entered  the  Earth's  atmosphere,  but  at  least  the 
Australian  tektites  with  few  exceptions  could  not  have  been  either  tumbling 
or  oscillating  with  large  amplitude  during  the  period  of  greatest  heating  or 
the  patterns  characteristic  of  ablation  on  bodies  with  fixed  orientation 
would  not  be  present  on  the  many  existing  specimens.  Therefore,  their  turn-, 
bling  motion  had  to  have  been  arrested  and  their  subsequent  oscillatory  motion 
reduced  to  quite  small  amplitudes  before  their  period  of  greatest  heating,  and 
it  is  the  consequences  of  this  requirement  that  can  be  studied  by  means  of  the 
present  analysis . 

Figures  of  equilibrium. -  In  order  to  apply  the  results  of  the  preceding 
sections,  a  value  of  initial  angular  velocity  must  be  assigned  to  the  body. 
Further,  the  body  shape  should  be  chosen  to  be  compatible  with  the  assigned 
value  of  angular  velocity,  for,  as  mentioned  above,  the  body  in  its  molten 
state  will  tend  to  assume  a  figure  consistent  with  its  turning  rate.  Let  it 
be  assumed  first  that  the  viscosity  of  the  body  in  its  molten  state  is  suffi¬ 
ciently  low  to  permit  the  stable  figure  to  be  attained  before  the  body  solid¬ 
ifies.4  The  form  attained  may  then  be  calculated.  In  effect,  the  problem  is 
analogous  to  the  famous  problem  in  cosmogony  of  determining  the  figure  of 
equilibrium  of  a  rotating  liquid  mass  (cf.,  e.g.,  ref.  12)  except  that  the 
force  tending  to  contract  the  mass  is  surface  tension  rather  than  gravitation 
(it  is  easy  to  show  that  for  bodies  the  size  of  tektites  the  mutual  attraction 

4It  is  probable  that  this  was  the  case  only  for  the  tektites  found  in  the 
part  of  the  strewnfield  covering  southeast  Australia  (cf.  ref.  11).  The  suc¬ 
ceeding  analysis  is  therefore  limited  to  that  particular  group  of  tektites. 
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of  the  bodies1  particles  is  negligible  compared  -with  the  surface  tension 
force).  The  definitive  analysis  is  due  to  Charrueau  (ref.  13) >  but  its  essen¬ 
tials  are  repeated  more  accessibly  in  volume  k  of  Appellfs  treatise  on  rational 
mechanics  (ref.  14) .  For  present  purposes,  it  suffices  to  note  the  following 
results: 

(1)  In  the  absence  of  gravitation,  the  figure  of  equilibrium  must  have 
cylindrical  symmetry;  that  is,  cross  sections  normal  to  the  axis  of  rotation 
are  circles. 

(2)  The  meridian  curve  for  the  figure  of  equilibrium  is  an  elliptic 
function  dependent  on  a  single  parameter  K2  having  the  form 

K2  =  Sto2a3/8f  (47) 


Sketch  (e) 


where 

5  liquid  density 

to  angular  velocity 

f  surface  tension 

a  radius  of  figure  in  equatorial 
plane  (see  sketch  (e)) 

(3)  For  K2  =  0,  that  is,  for 
zero  turning  rate,  the  figure  of  equi¬ 
librium  is  a  sphere.  Let  its  mass  be 
m  and  its  radius  R.  A  body  of  the 
same  mass  with  a  small  turning  rate 
will  flatten  at  its  poles.  For  K2«l; 
it  is  essentially  an  ellipsoid  of 
revolution.  The  meridian  curve  z(r) 
(see  sketch  (e))  is  approximately 

z  »  (1  -  K2)  Ja2  -  r2  ;  K2  «  1 

(48) 


where,  for  the  same  mass  m,  a  and  b,  the  semiaxes  of  the  ellipse,  are 
related  to  R  by 


(49) 
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(5)  In  the  absence  of  an  external  pressure,  there  are  no  figures  of 
equilibrium  for  K2  >  1.  Presumably  in  this  case  the  body  "will  decelerate 
through  a  progression  of  elongating  pear-shaped  figures  culminating  in 
separation  at  the  weakest  section. 


In  application  to  tektites,  it  remains  to  assign  representative  numeri¬ 
cal  values  to  the  physical  properties  appearing  in  the  parameter  K2.  These 
may  be  obtained  from  results  given  in  reference  4.  The  density  5  of  tek- 
tite  glass  is  given  as  2.4  gm/cm3.  The  surface  tension  f  for  glass  of  tek- 
tite  composition  is  reported  as  3^0  dynes/cm.  A  single  representative  mass 
will  be  assumed  and  will  be  taken  to  be  that  for  a  sphere  of  1-cm  radius. 

With  these  numbers,  the  limiting  values  of  co  and  angular  momentum  for  which 
cylindrical  symmetry  is  sought  (K2  =  l)  are,  respectively,  the  order  of 
25  radians /sec  and  l6j  dyne-cm-sec. 


Probability  of  a  tumbling  entry.-  It  has  been  assumed  that  the  body 
attains  a  figure  of  equilibrium  in  space  before  it  solidifies.  As  a  conse¬ 
quence,  it  must  rotate  at  a  constant  angular  velocity  about  an  axis  through 
its  center  of  gravity,  and  that  axis  must  maintain  a  fixed  attitude  with 
respect  to  space-fixed  coordinates  (ref.  14) .  Hence,  the  angular  velocity 
vector  also  maintains  a  fixed  inclination  with  respect  to  the  velocity  vector 
and  the  body  enters  the  atmosphere  in  this  condition.  An  inclination  of 


entry  as  defined  here,  whereas  an  incli- 

Con- 


exactly  90°  corresponds  to  a  tumbling 

nation  of  0°  corresponds  to  a  rolling  motion  around  the  velocity  vector 

sider  first  the  consequences  of  a 
purely  tumbling  entry  for  a  figure  of 
equilibrium,  that  is,  for  a  body  having 
cylindrical  symmetry  about  the  axis  of 
rotation,  the  latter  being  inclined  90° 
to  the  stream.  As  shown  on  sketch  (f), 

I _  the  aerodynamic  force  on  the  body 

v  always  remains  parallel  to  the  velocity 
vector  and  passes  through  the  bodyTs 
center  of  volume.  Barring  the  exist¬ 
ence  of  a  significant  inhomogeneity  or 
asymmetry,  the  center  of  gravity  is 
coincident  with  the  center  of  volume 
and  the  aerodynamic  restoring  moment 
about  the  center  of  gravity  is  identi¬ 
cally  zero.  Hence,  the  body  in  this  mode  is  in  a  state  of  unstable  equilib¬ 
rium;  a  disturbing  force  whose  moment  vector  is  even  slightly  misalined  with 
the  angular  velocity  vector  will  cause  the  body  to  depart  from  this  mode  and 
seek  a  new  state  of  motion  in  which  it  is  stable.  Now  consider  a  case  in 
which  the  axis  of  rotation  is  inclined  to  the  velocity  vector  at  an  angle  less 
than  90°.  Sketch  (g)  shows  a  view  of  the  body  in  the  plane  containing  both 
the  axis  of  rotation  and  the  velocity  vector.  As  the  body  has  cylindrical 
symmetry  about  the  axis  of  rotation,  its  projection  in  this  plane  is  invariant 
with  time;  the  profile  is  approximately  elliptical  as  a  result  of  the  flatten¬ 
ing  the  body  has  undergone  while  turning  in  the  liquid  state.  The  aerodynamic 
force  on  the  body  lies  wholly  and  continually  in  the  plane;  it  is  approxi¬ 
mately  alined  with  the  stream,  but  passes  through  a  point  representative  of 
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the  center  of  curvature  of  the  forward 
surface,  a  point  which  lies  behind  the 
center  of  gravity.  It  is  clear  that 
in  this  case  an  aerodynamic  restoring 
moment  is  developed  about  the  center  of 
gravity  that  will  drive  the  angle  of 
attack  a  toward  zero.  These  consid¬ 
erations  make  it  evident  that  the  prob¬ 
ability  that  a  figure  of  equilibrium 
will  undergo  a  tumbling  entry  is  very 
remote.  The  stable  state  is  with  the 
cylindrical  surface  broadside  to  the 
stream,  and  for  all  initial  inclinations 
other  than  exactly  90°,  the  body  will 
tend  to  this  state. 


Nonplanar  motions. -  With  tumbling 
effectively  eliminated  as  the  normal  Sketch  (g) 

mode  of  entry  for  bodies  which  are  fig¬ 
ures  of  equilibrium,  consideration  must 
be  given  to  nonplanar  motions  in  which 

the  angular  velocity  vector  is  initially  fixed  in  space,  inclined  to  the 
velocity  vector  at  some  angle  less  than  90°  •  It  Is  possible,  however,  to 
relate  this  study  in  a  simple  fashion  to  results  already  obtained  here. 


In  reference  15,  Leon  treated 
the  problem  of  a  spinning  body  entering 
the  atmosphere,  but  under  the  restriction 
that  the  initial  inclination  of  the 
angular  velocity  vector  from  the  velocity 
vector  be  small.  As  initial  inclinations 
up  to  90°  are  of  interest  here,  the 
results  of  reference  15  are  of  only 
limited  applicability  in  the  present  con¬ 
text.  A  means  of  removing  this  limita¬ 
tion  has  been  found;  however,  as  it  is 
suggested  by  the  form  of  the  results  pre¬ 
sented  in  reference  15,  these  will  be 
reviewed  briefly  before  their  extension 
is  presented.  A  sketch  of  the  relevant 
coordinates,  adopted  from  reference  15, 
is  shown  in  sketch  (h) .  The  most  impor¬ 
tant  angle  is  a ,  the  inclination  of  the 
body  axis  z  about  which  the  body 
rotates,  to  the  velocity  vector  along 
Z.  As  the  body  has  cylindrical  symmetry 
about  the  z  axis,  the  aerodynamic  force 
on  the  body  is  a  function  only  of  a  and 
lies  in  the  cr  plane.  The  aerodynamic 
moment  about  the  center  of  gravity  there¬ 
fore  lies  on  an  axis  normal  to  the  o 
plane,  or  along  the  line  of  nodes.  The 


z 


Z  Flight  path  direction 
z  Body  axis  of  rotation 
a  Resultant  angle  of  attack 

Sketch  (h) 
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angle  a  in  effect  defines  the  move¬ 
ment  of  the  nose  of  the  body  about  the 
velocity  vector.  In  space  and  on  entry 
into  the  atmosphere  cr  is  fixed  in  k 
value  and  orientation,  but  as  the  aero¬ 
dynamic  moment  grows,  the  a  plane 
will  begin  to  rotate  about  Z  causing 
the  body  nose  to  move  in  a  diminishing 
spiral.  A  typical  a  history  is 
shown  on  sketch  (i)  both  as  it  appears 
to  an  observer  on  the  Z  axis  and  as 
a  function  of  time.  A  very  useful 
relation  is  found  in  reference  15, 
namely,  the  ratio  between  the  asymp¬ 
totic  behavior  of  .the  envelope  of  cr, 
(a)env>  (cf.  sketch  (i))  to  the  asymptotic  behavior  of  the  envelope  of  a  for 
a  planar,  nonspinning,  nontumbling  entry.  For  the  same  entry  conditions 
(cfi  =  cxq,  a±  =  aq  =0),  the  ratio  is  simply 


(°')env 
(a)  env 

where 


Thus,  it  is  indicated  that  the  ratio  of  the  envelopes  is  a  constant  and  is 
increased  by  a  factor  dependent  only  on  a  single  parameter  v. 

The  simple  form  of  equation  (51)  suggested  the  possibility  that  a  similar 
form  would  result  even  with  the  removal  of  the  limitation  to  small  values  of 
oq.  This  has  been  found  to  be  the  case.  When  the  aerodynamic  restoring- 
moment  coefficient  can  be  expressed  as  sin  cr,  the  extension  of  equa¬ 

tion  (51)  to  the  case  of  arbitrarily  large  values  of  a±  is  of  the  form 

(tan  cr/2)env  /  rt7/2  G(a±,v)  .  -  2  ^  . 

(tan  oo/2)env  J  tanh  jtv/2  ^  ) 


7tV 


/2 


tanh  (jtv/2) 


oo 


(51) 


Y 


In  the  present  study,  only  small  to  moderate  values  of  v  need  be  considered, 
ancl  in  this  case  the  ratio  involving  G  is  essentially  unity.  Then,  as 
before,  the  ratio  of  the  envelopes  is  increased  by  a  simple  factor  dependent 
only  on  v.  As  the  planar  solution  is  already  available  (eq.  (30)),  equa¬ 
tion  (52)  provides  a  simple  means  of  studying  the  angle -of -attack  behavior  of 
rotating  bodies  having  cylindrical  symmetry  about  the  axis  of  rotation.  In 
effect,  one  need  only  study  the  simplest  of  planar  problems,  nontumbling 
entry  from  without  the  atmosphere;  multiplying  the  result  by  the  factor  in 
equation  (52)  then  gives  the  desired  property  of  the  rotating  body. 


In  this  light ,  consider  the 
nontumbling  planar  entry  of  the  figure 
of  equilibrium  with  its  cylindrical  sur¬ 
face  facing  the  stream.  On  sketch  (j), 
a  figure  of  equilibrium  is  shown  super¬ 
imposed  on  a  sphere  of  the  same  mass  so 
that  their  front  faces  are  alined  as 
nearly  as  possible.  If  the  front  face 
of  the  figure  of  equilibrium  could  be 
alined  perfectly  with  that  of  the  sphere, 
it  is  clear  that  the  moment  about  the 
center  of  gravity  of  the  figure  of  equi¬ 
librium  would  be  of  the  form 


=  -Cp 


sin  a 


Sketch  (j) 


It  is  indicated,  therefore,  that  for 
small  amounts  of  oblateness  equation  (53) 

should  adequately  represent  the  pitching  moment  for  the  figure  of  equilibrium. 
On  the  other  hand,  because  of  the  fore-and-aft  symmetry  of  the  figure,  the 
pitching  moment  must  be  zero  at  a  =  rt/2,  so  that  equation  (53)  will  not 
apply  for  values  of  a  in  the  immediate  vicinity  of  jr/2.  With  this  range 
excluded,  equation  (53)  is  applicable  to  the  figure  of  equilibrium  having 
small  oblateness.  Hence,  with  Cm  of  the  required  form,  equation  (30)  is 
applicable  and  is  rewritten  here  for  convenience 


(tan  e)  ( 
tan  e 


=  f^G( 

sjj XX 


where  e  =  -a/2.  It  is  convenient  to  rewrite  equation  (54)  in  terms  of -the, 
ratio  q/q^,  where  q^  is  the  dynamic  pressure  at  maximum  heating.  We  have 
from  equation  (32) 


X2  =  K2 


. W  1 


From  equations  (35)  and  (4-1) 


<imh  =  \  e'1/3yi2P sin  n  (cy) 

so  that  x2  becomes,  when  equations  (10),  (53)>  (55)  >  and  (50  are  used, 


_  £  -1/3 

"3 


P  sin  y 


(R  -  ,b) 


As  only  small  amounts  of  oblateness  are  admissible,  the  figure  of  equilibrium 
is  essentially  an  ellipsoid  of  revolution  for  which  the  expressions  for  m  and 
b  previously  presented  (eqs.  (47)  through  (50))  may  be  used.  Retaining  terms 
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only  to  the  first  order  in  K2  gives 


(R  '  b)  f  *  f  f  *  K2  «  1  (58) 

which  relates  the  geometric  properties  of  the  body  to  the  rotation  parameter 
K2.  The  envelope  expression  (eq.  (54))  becomes 

(r 6>env  -  Fo(^) - i -  (59> 

tan  €-p  \/  ic  p'  .1/4  i/4 

fk  e-i/3 

\9  PR  sin  7j_/ 


Finally,  the  envelope  expression  for  the  initially  inclined  rotating  body, 
(tan  cr/2)env,  is  obtained  simply  from  the  product  of  equation  (59)  and  the 
factor  given  in  equation  (52).  The  result  is 


(tan  cr/ 2)  0ny 
tan  c p 


K  tanh 


JtV 


G(ep) 


PR  sin 


1/4 


(60) 


Equation  (60)  reveals  an  interesting  result.  It  will  be  noted  that  both  K 
(eq.  (47))  and  v  (eq.  (51))  are  linear  functions  of  w.  As  tanh(rtv/2)  is 

essentially  unity  for  all  values  of  V  greater  than  about  1.5*  equation  (60) 

indicates  that  the  effects  on  the  envelope  due  separately  to  K 2  and  v  cancel 
each  other.  On  the  one  hand,  increasing  K2  corresponds  to  increasing  the 
turning  rate  of  the  liquid  body,  giving  it  a  greater  oblateness,  and  this  in 

turn  leads  to  a  larger  restoring  moment  which  makes  the  envelope  subside  more 

rapidly.  On  the  other  hand,  the  inertial  effect  of  the  larger  turning  rate  is 
to  widen  the  envelope.  These  two  opposing  effects  cancel.  The  result  is 
that,  regardless  of  the  turning  rate,  all  bodies  formed  by  turning,  having  the 
same  mass  and  the  same  initial  inclination  cr  ,  will  have  identical  angle- of- 
attack  envelopes.  Although  this  result  has  been  obtained  under  the  restric¬ 
tion  of  small  turning  rates,  numerical  solutions  for  larger  rates  indicate 
that  it  remains  essentially  true  for  all  rates  up  to  the  limiting  value 
(K2  =  1)  for  which  a  figure  of  equilibrium  exists. 


Consider  the  application  of  this  result  to  the  southeast  Australian 
tektites.  Of  all  the  round-form  tektites  that  have  been  discovered  in  this 
region,  none  have  been  reported,  to  the  authors T  knowledge,  to  have  any  more 
than  a  small  amount  of  oblateness;  that  is,  in  reference  to  figure  9>  none 
apparently  have  primary  figures  corresponding  to  values  of  K2  greater  than 
perhaps  0.2.  This  means  that  if  they  had  been  formed  by  turning  in  the  liquid 
state,  their  turning  rates  had  to  have  been  small.  The  analysis  above  leads 
one  to  expect  that,  as  the  tektites  having  small  turning  rates  evidently  were 
able  to  survive  their  passage  through  the  atmosphere,  tektites  having  forms 
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consistent  with  considerably  larger  rates  should  have  been  able  to  survive  as 
well  since  their  angle -of -attack  envelopes  should  have  been  essentially  the 
same.  That  such  tektites  apparently  have  not  been  found  in  southeast 
Australia  therefore  cannot  be  attributed  to  their  having  failed  to  survive 
their  flight  through  the  atmosphere,  but  must  be  charged  to  the  mechanism  by 
which  the  tektites  were  formed.  The  evidence  is,  then,  that  this  mechanism 
was  capable  of  imparting  to  the  tektites  only  a  very  limited  amount  of 
angular  momentum. 

Effect  of  ablation.-  The  analysis  of  the  preceding  section  should  ade¬ 
quately  describe  the  behavior  of  the  round -form  tektites  over  the  initial  por¬ 
tion  of  their  entry  into  the  atmosphere.  As  soon  as  ablation  begins,  however, 
the  bodies  undergo  a  marked  change  in  geometry,  and  this  must  be  taken  into 
account  in  the  analysis  of  their  subsequent  behavior.  This  may  be  done  quali¬ 
tatively  within  the  framework  of  the  present  theory  in  view  of  the  following: 
First,  as  in  the  preceding  section,  one  need  only  study  the  case  of  planar 
motion  and  then  multiply  the  result  by  the  factor  given  in  equation  (52)  to 
obtain  the  behavior  of  the  rotating  body.  Second,  in  consideration  of  the 
planar  motion,  it  may  be  said  that  the  previous  analysis  will  hold  over  the 
range  of  dynamic  pressure  from  essentially  zero  to  a  value  q^  at  which  abla¬ 
tion  begins.  Results  of  the  ablation  studies  of  reference  4  indicate  that  the 
ratio  q^/q^ih  varies  somewhat  with  entry  conditions  and  tektite  size.  For 
entry  conditions  compatible  with  a  lunar  origin  and  for  tektites  the  size  of 
those  found  in  southeast  Australia,  results  from  reference  4  indicate  a  value 
for  the  ratio  of  approximately  l/20.  For  convenience,  let  it  be  assumed  that 
a  peak  in  the  oscillatory  motion  occurs  in  the  vicinity  of  this  point..  The 
corresponding  values  of  €p  and  x(ep)  then  may  be  determined  from  the  preced¬ 
ing  analysis  and  will  serve  as  the  initial  conditions  for  the  ensuing  motion. 

As  ablation  begins,  the  surface 
facing  the  stream  will  begin  to  recede, 
increasing  its  radius  of  curvature  (cf. 
fig.  19,  ref.  4).  It  is  reasonable  to 
assume  that,  though  the  angle  of  attack 
continually  changes,  the  forward  surface 
will  continue  to  present  an  essentially 
spheroidal  face  to  the  stream.  Hence, 
approximately,  the  aerodynamic  force  on 
the  body  will  remain  alined  with  the 
stream  direction  and  pass  through  the 
center  of  curvature  of  the  forward  sur¬ 
face.  This  is  illustrated  on  sketch  (k) . 

As  the  forward  surface  continues  to 
recede,  its  center  of  curvature  moves 
rearward.  The  center  of  gravity  also 
moves  rearward,  but  less  rapidly  than  the 
center  of  curvature,  so  that  the  static 
margin  £  increases  with  time.  The 
result  is  that  the  aerodynamic -restoring 
moment  about  the  center  of  gravity  retains 
the  form 

Cm  =  -Cp|  sin  a 


Sketch  (k) 


(61) 
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but  now  Cd  and  especially  |  are  increasing  functions  of  time.  The  equa¬ 
tion  of  motion,  equation  (4),  also  retains  the  same  form 


a  -  ^  Cp|q(t)  sin  a  =  0  (62) 

where  now  I  also  is  a  function  of  time.  However,  for  the  spheroidal  body 
the  changes  in  I  and  Cp  are  small  enough  to  be  neglected,  or,  in  any  case,  I 
and  CD  may  be  combined  with  £  to  give  an  "effective"  value  of  that  param¬ 
eter.  Finally,  since  both  £  and  q(t)  increase  with  time,  it  may  be  assumed 
that,  at  least  qualitatively,  |  increases  as  some  power  of  q,  that  is, 


I 


(63) 


where  £a  and  qa  are,  respectively,  the  static  margin  and  dynamic  pressure  at 
the  beginning  of  ablation.  On  insertion  of  the  exponential  approximation  for 
q  (eq.  (8)),  the  equation  of  motion  becomes 

a  -  q„  Cp£  eS(1+A)t  sin  a  =  0  (64) 

d  ct 


which  retains  the  form  for  the  Painleve  transcendent.  In  effect,  the  increas¬ 
ing  static  margin  causes  the  body  to  behave  as  though  it  were  a  nonablating 
body  passing  through  an  atmosphere  with  a  larger  density  gradient  than  that 
actually  existing. 


With  the  form  of  the  Painleve  equation  retained,  all  of  the  results 
obtained  previously  for  the  planar  oscillatory  motion  may  be  made  applicable 
merely  by  a  change  in  notation.  Thus,  let 

s  =  s(l  +  A)  1 


A  2  ^  u 

—  cDla 


X 


—  st/2 
=  Ke 


(65) 


J 

The  asymptotic  behavior  of  the  envelope  curve  follows  from  equation'  (28) 


(tan  e)eny  C(k)g(£P,k)  . 

tan  ep  ”  s/f  . 

It  is  useful  to  cast  this  result  in  terms  of  the  dynamic  pressure  ratio  q/qmh. 
where,  as  before,  q^  is  the  dynamic  pressure  at  maximum  heating.  The  result 
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is 


<tan  e>env 
tan  Gp 


C(K)G(ep,K) 


(l+A)/4/  ^  \(  l+A)/4 
V^nh/ 


(67) 


The  parameter  k  may  be  related  to  the  original  coordinate  x  at  the  begin¬ 
ning  of  ablation  through  the  expression 


K  = 


x  (6p) 

1  +  A 


(68) 


Finally,  the  complete  envelope  curve  for  the  planar  motion,  obtained  from 
equations  (59)  and  (67)  is  multiplied  by  the  factor  J  (rtv/2)/ [tanh(jr v/2)]  to 
give  the  desired  result  (tan  cr/2)env.  Results  obtained  from  the  above  analy 
sis  are  shown  on  figure  10  for  several  values  of  the  ablation  parameter  A. 


The  body  has  the  mass  of  a  sphere  of 
1-cm  radius  and  a  turning  rate  of 
1  radian/sec.  Entry  conditions  are 
those  estimated  in  reference  4  to  be 
compatible  with  a  lunar  origin: 

Vi  =  11.2  km/ sec 
P  =  0.1396  km"1 
7±  =  20° 


The  main  results  are  considered  to  be 
those  shown  on  figure  10(a)  for  an 
initial  inclination  a±  of  40°  since 
if  the  initial  inclination  was  truly 
arbitrary,  the  greatest  number  of  tek- 
tites  should  have  entered  the  atmos¬ 
phere  with  inclinations  near  the  mean  -g 
between  0°  and  90°.  Figure  19  of  refer-  > 
ence  4  was  used  to  establish  a  represent-  b 
ative  value  of  A  for  tektites  undergoing 
ablation  under  nonosc dilatory  conditions: 
Over  the  important  initial  period  of. 
ablation,  during  which  time  the  radius  of 
curvature  of  the  front  face  increases 
very  rapidly,  A  was  found  to  be  of  the 
order  of  4  or  5*  For  the  initially 
inclined  rotating  tektite,  of  course  a 
smaller  value  must  be  expected,  but  it  is 
seen  on  figure  10(a)  that'  even  a  fraction 
of  the  full  amount  is  sufficient  to  reduce 


0  .2  .4  .6  .8  1.0 


q/q*h 

(b)  a±  =  8o° 

Figure  10.-  Envelopes  of  resultant  angles  of 
attack  of  a  rotating  spheroidal  tektite 
for  a  range  of  values  of  the  ablation  . 
parameter  A;  turning  rate 
co  =  1  radian/sec. 
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a  to  negligible  proportions  well  before  maximum  heating.  The  results  lead 
to  the  conclusion  that  entry  conditions  compatible  with  a  lunar  origin  permit 
sufficient  time  for  the  greater  part  of  the  round-form  Australian  tektites  to 
aline  themselves  with  the  stream  well  before  maximum  heating,  in  compliance 
with  the  physical  evidence  of  their  having  undergone  ablation  under  nonoscil- 
latory  conditions  over  the  high  heating  portion  of  their  trajectories.  Shown 
on  figure  10(b)  are  results  for  an  initial  inclination  of  80°.  It  is  noted 
that,  if  the  same  rate  of  ablation  is  allowed  (A  =  3)  as  that  which  brought 
the  body  to  negligible  amplitudes  in  the  case  of  a  k0°  initial  inclination, 
here  the  body  still  retains  a  residual  amplitude  as  maximum  heating  is 
approached,  but  of  no  more  than  a  few  degrees.  Then  in  this  case  also,  the 
final  stage  of  ablation  occurs  with  the  body  in  an  essentially  fixed  attitude 
with  respect  to  the  stream,  so  that  the  ring  waves  and  coiled  flanges  charac¬ 
teristic  of  this  stage  should  make  their  appearance  on  this  body  as  well.  On 
the  other  hand,  since  over  the  initial  portion  of  the  ablative  process  the 
inclination  of  the  nose  of  the  body  was  considerably  greater  than  in  the  first 
case,  the  nose  will  have  received  considerably  less  heating  than  in  the  first 
case,  and  hence  will  have  receded  a  lesser  amount.  The  initial  inclination 
must  be  considered  to  be  arbitrary  and  therefore,  for  initially  identical 
bodies,  any  amount  of  ablation  is  possible  between  the  maximum,  when  the  ini¬ 
tial  inclination  is  zero,  and  that  for  initial  inclinations  near  90°.  It  is 
believed  that  this  explains  how  it  can  happen  that  initially  identical  tek¬ 
tites  following  identical  trajectories  may  nevertheless  show  quite  different 
depths  of  recession  of  their  forward  surface  (cf.  ref.  4). 

Nonspheroidal  shapes. -  While  the  preceding  analysis  would  appear  to  apply 
to  the  majority  of  the  tektites  found  in  southeast  Australia  (i.e.,  those  of 
this  region  whose  shapes  are  essentially  spheroidal),  it  does  not  account  for 
the  remaining  shapes,  in  particular  the  ellipsoids  and  "dumbbells."  Figure  11 
Illustrates  the  three  shapes,  the  original  unablated  forms  being  evident  in  the 
rear  view.  It  is  interesting  to  speculate  on  possible  mechanisms  that  might 
account  for  the  latter  two  forms.  The  fact  that  oblate  spheroids,  ellipsoids, 
and  dumbbells  make  their  appearance  also  in  the  study  of  rotating  liquids 
when  the  contractive  force  is  gravitation  (ref.  12)  suggests  at  first  glance 
that  the  mechanism  for  the  formation  of  tektites  might  have  been  of  a  similar 
nature.  The  investigation  of  this  possibility  was  in  fact  what  motivated  the 
present  research.  In  the  gravitational  case,  the  spheroids  and  ellipsoids 
form  a  progression  of  figures  of  equilibrium  with  increasing  angular  momentum, 
while  the  dumbbell  or  pear-shaped  figure  makes  its  appearance  for  values  of 
angular  momentum  beyond  that  of  the  last  stable  ellipsoid  (whether  the  pear- 
shaped  figure  is  itself  stable  is  a  delicate  and  still  controversial  question). 
Unfortunately,  the  case  for  the  existence  of  a  parallel  phenomenon  for  the 
tektites  is  untenable.  When  the  contractive  force  is  surface  tension  the 
figure  of  equilibrium  must  have  cylindrical  symmetry  about  the  axis  of  rota¬ 
tion  (ref.  1^),  which  rules  out  both  the  ellipsoid  and  the  dumbbell. 

Two  alternative  phenomena  may  be  envisaged:  First,  these  bodies  may  have 
had  angular  momenta  large  enough  to  exceed  the  limiting  value  for  which  a 
stable  figure  exists.  They  would  then  have  sought  to  rid  themselves  of  part 
of  their  angular  momentum,  and  this  process  could  conceivably  have  led  to  the 
formation  of  elongated  ellipsoids  and  dumbbells,  these  solidifying  before  they 


36 


could  break  apart.  Let  us  confine  attention  only  to  specimens  of  a  common 
initial  mass,  and  let  the  mass  be  that  of  a  sphere  of  1-cm  radius,  a  charac¬ 
teristic  shared  by  a  large  number  of  the  specimens  found  in  southeast 
Australia.  Then  according  to  the  above  hypothesis,  the  elongates  must  have 
had  angular  momenta  in  excess  of  approximately  170  dyne-cm-sec.  On  the  other 
hand,  spheroids  of  the  same  mass  found  in  southeast  Australia  could  have  had 
angular  momenta  only  of  the  order  of  0-60  dyne-cm-sec  on  the  presumption 
already  noted  that  their  values  of  K2  were  no  larger  than  perhaps  0.2.  Then 
there  is  a  range  of  angular  momenta  from  60  to  I70  dyne-cm-sec  in  which  fig¬ 
ures  of  equilibrium  consistent  with  values  in  this  range  apparently  have  not 
been  found  in  southeast  Australia.  As  has  been  discussed,  the  dynamic  behav¬ 
ior  of  these  bodies  should  have  been  similar  to  that  of  surviving  bodies  hav¬ 
ing  smaller  values  of  angular  momentum,  and  hence,  there  is  no  evident  reason 
why  they  should  have  failed  to  appear.  That  is,  if  the  elongates  had  been 
formed  as  the  result  of  having  an  excess  of  angular  momentum,  then  flattened 
spheroids  should  have  been  formed  as  well  and  should  have  been  found  in  south¬ 
east  Australia.  This  inconsistency  leads  one  to  doubt  that  the  elongates 
could  have  had  such  large  values  of  angular  momentum  and  therefore  to  reject 
the  idea  that  they  could  have  been  formed  in  the  way  envisaged  above. 

With  rotation  rejected  as  the  mechanism  by  which  the  elongates  acquired 
their  form,  a  second  alternative  is  the  break-up  of  a  slowly  turning  jet  of 

liquid  glass.  As  is  well  known 
(ref.  l6),  a  jet  of  fluid  is  unstable 
and  will  exhibit  ever-increasing  har¬ 
monic  undulations  along  its  length.  As 
shown  in  the  upper  section  of 
sketch  (z),  a  freed  drop  would  tend  to 
exhibit  an  ellipsoidal  form,  whereas, 
as  in  the  lower  section  of  the  sketch, 
the  drop  would  tend  to  the  dumbbell 
form.  Surface  tension  would  then  cause  all  forms  to  tend  to  the  spheroidal. 
The  appearance  of  all  three  forms  on  the  surface  of  the  Earth  may  be  attrib¬ 
uted  to  varying  rates  of  solidification,  some  bodies  having  solidified  before 
attaining  their  final  forms . 


o 


o-c=o 


Sketch  (z) 


CONCLUDING-  REMARKS 


The  tumbling  motion  of  aerodynamically  stable  bodies  entering  planetary 
atmospheres  has  been  analyzed  considering  that  the  tumbling  motion,  its  arrest, 
and  the  subsequent  oscillatory  motion  are  governed  by  the  differential  equation 
for  the  fifth  Painleve  transcendent.  A  study  of  the  asymptotic  behavior  of  the 
transcendent  enabled  the  functional  relationship  between  the  envelope  of  oscil¬ 
latory  motion  and  all  the  significant  body  and  planetary  properties  to  be  demon¬ 
strated  in  a  concise  expression.  Results  were  applied  to  the  study  of  vehicles 
intended  for  use  as  planetary  probes  and  it  was  shown  how  rapid  estimates  could 
be  made  of  their  probable  amplitudes  of  oscillation  in  relation  to  aerodynamic 
heating  and  loads.  The  theory  was  also  applied  to  a  study  of  the  motions  of  the 
southeast  Australian  tektites.  It  was  concluded  that  with  entry  conditions 
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compatible  with  those  for  a  lunar  origin,  a  considerable  range  of  initial 
turning  rates  and  initial  inclinations  is  admissible  for  the  spheroidal  tek- 
tites  within  which  their  amplitudes  of  oscillatory  motions  are  reduced  to. 
negligible  proportions  before  maximum  heating.  This  complies  with  the  evi¬ 
dence  of  their  having  undergone  ablation  at  essentially  fixed  attitudes,  oyer 
the  high  heating  portion  of  their  trajectories.  The  admission  of  an  initial 
inclination  from  the  direction  of  flight  of  the  axis  about  which  the  body 
turns  leads  to  an  explanation  of  how  it  can  happen  that  initially  identical 
tektites  following  identical  trajectories  may  show  significantly  different 
depths  of  ablative  recession  of  their  forward  surfaces.  It  is  suggested  that 
the  mechanism  by  which  the  southeast  Australian  tektites  acquired  their  forms 
was  the  break-up  of  a  slowly  turning  jet  of  liquid  glass. 
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"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof ” 

— National  Aeronautics  and  Space  Act  of  1958 
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